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^ , Abstract 

In — > i/ M oscillations, we find that the effect of the CP phase S becomes large in the region E < 2 
GeV and L > 2000 km. In this region, the change of the probability in this channel reaches about 0.4 
due to the CP phase effect beyond our expectation in the case of large 1-3 mixing angle. Furthermore, 
the CP phase effect have almost same sign over the region E > 0.5 GeV so that one may find the 
signal of CP violation by measuring the total rate only. As an example, we use an experimental setup 
and demonstrate that the allowed region is limited to one by combined analysis of u e and f M events 
^ , although there remain three allowed regions by the analysis of v e events alone. 
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^ ! 1 Introduction 

\ The measurement of the leptonic CP phase is one of the most important aims in elementary particle 

■ physics. A lot of investigations for this possibility have been performed, see [TJ [21 El HI El El [3 [S] and the 

references therein. In most of these investigations, v e — > or — > v e oscillations (so-called golden channel 
[2]) are used for the measurement of 6 and z/ e — > v r oscillations (silver channel [5]) are also used partly. 
However, there are also — > oscillations as another channel observed in long baseline experiments. 
The measurement of 5 in this channel have been considered to be difficult because the probability 
depends almost only on the cos<5 term and the CP dependence disappears in the case of 5 = 90° and 
270°. However, some recent papers discuss the possibility of solving the parameter ambiguity by using this 
channel [9l [10] . In these analysis, the CP dependence becomes small and the confirmation of 5 seems to be 
difficult in actual experiments. However, if the CP dependence becomes large in this channel, this analysis 
is very useful, interesting and cost-effective because the probability of this channel can be measured at the 
same time as — » v e probability without any extra apparatus and we can obtain another information on 
the CP phase in superbeam experiments. In quark sector, the CP violation was discovered in K physics, 
but the determination of the CP phase cannot be performed by the increase of the data in K physics and 
needed the B-factory experiments [TT]. The measurement of the CP phase in two different channels is 
considered to be very important also for the lepton sector. We obtain the information with a different 
nature by using different channels. Furthermore, this opens the window for a verification of the unitarity 
in three generation and exploration of new physics. 

In the previous paper, we introduced a new index of the leptonic CP phase dependence lev and 
investigated how the information of 5 can be obtained through the channel of — > v e oscillation [12] . In 
this letter, we use Icp for the channel of — > oscillation and explore the region in E-L plane where the 
CP dependence becomes large. We consider one baseline in this region and investigate the contribution 
of the CP phase to — > oscillations in detail. Finally, we simulate both and v e events by using an 
experimental setup and compare the analysis of v e events alone with the combined analysis of and v e 
events. 



2 Large CP Dependence in — > Oscillations 

:r [T3] as 

H = 23 TH'^Oj 3 (1) 



At first, we write the Hamiltonian in matter 13 as 
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by factoring out 9 23 and S, where 23 is the rotation matrix in the 2-3 generations and T is the phase 
matrix defined by T — diag(l, 1, e iS ), The reduced Hamiltonian H' is given by 

H' = Oi 3 Oi 2 diag(0, A21, A 3 i)Of 2 Of 3 + diag(a, 0, 0), (2) 

where A, tj = Amf J /(2E) = (to 2 - mf)/(2E), a = y/2G F N e ~ 7.56 x 1CT 5 • pY e , G F is the Fermi constant, 
N e is the electron number density, Y e is the electron fraction, E is neutrino energy and rrii is the mass 
of Vi. The oscillation probability for v u — > v a is proportional to the cos S and cos2<5 in constant matter 
profile [14] and can be expressed as 

P m = A cos S + C + D cos 25. (3) 
Here, A, C and D are determined by parameters other than 6 and are calculated by 

A = 4Rc[(^^3 + S' TTS 2 23 yS^}c 23 s 23: (4) 

c = \s'j 2 4 3 + \s> T \ 2 s i 3 (5) 

+2(\S' nT \ 2 +Re[S^S' TT ])c 2 23S 2 23 , ( 6 ) 

D = 2\S'^cl 3 sl 3 , (7) 

where S' a p = [exp(— iH'L)] a ^. In the case of symmetric density profile and 6> 23 = 45°, the coefficient A 
coincides with A eil , which is the coefficient of cos S in the — > v e oscillation probability, except for the 
sign [TS]. Namely, A is simplified as 

A=-A e , = -Re[S'; e S' Te }. (8) 

If we notice the small parameters, a = A 2 i/A 31 and sin^ = s 13 and count the order of coefficients, we 
obtain A — 0(as 13 ), C — 0(1) and D = 0(a 2 s\ 3 ). Therefore, this leads to the relation D <C A <C C and 
we can safely neglect the term including D [T5] . 

Next, we introduce Icp following |12j . Suppose that P m ax and P m i n as the maximal and minimal values 
when 5 changes from 0° to 360°. Then, Icp is defined by 

T Pmax Pmin / n \ 

/cp = p Tp— (9) 

J max ^ ± mm 

also in — > v u oscillation. As the coefficient D can be neglected, P MM is approximated by P MAI ~ A cos 5+C. 
Pmax and P m in correspond to the two end points. Namely, the maximal and minimal values are given by 
Pan — A + C or — A + C corresponding to S = 0° or 180°. So, we obtain 

^cp^- (10) 

Next, let us present the region with large Icp in E-L plane. In order to exclude the region with large 
7qp due to the small denominator, we also show the region with large \A\, which is the numerator of /cp- 
In fig.l, we assume sin 2 26>i 3 = 0.16, which is near the upper bound of the CHOOZ experiment |16j . and 
can be measured in the next generation reactor experiments [TTJ [TS] . If #i 3 will be determined by this 
method without ambiguity, we can take a large step for the measurement of 5 [19]. We use the best-fit 
values Am^ = 7.9 x l(T 5 eV 2 , sin 2 12 = 0.31, Amf^ = 2.2 x 10~ 3 eV 2 and sin 2 23 = 0.50 [5D] for other 
parameters. We calculate by using the matter density p = 3.3g/cm 3 . In fig.l, black color shows the region 
with Icp > 30% and |A| > 20% respectively in left and right figures. Comparing the two figures, we can 
see that the region with large Iqp does not overlap that with large |^4| above 2 GeV. This means that the 
region of large Iqp above 2 GeV just indicates the region where the denominator is small. On the other 
hand, both regions overlap over wide region below 2 GeV. We found from fig.l that the CP dependence 
in — * oscillations becomes large in the region E < 2 GeV and L > 2000 km. We choose a baseline 
of L = 3000 km as an example in this region and investigate the behavior of A and P MAt in detail. See ref. 
[2"T] about other baselines. 
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Fig.l. Region with large Iqp (Left) and large \A\ (Right). The region with black color has a value larger 
than 30% for Icp and 20% for |j4|. Left and right panels show the region for Iqp and 

In fig. 2, we calculate the energy dependence of A and P MM at L = 3000 km by using the exact expressions 
([5]) and ([3]). P MM is plotted in the energy region E — 0.4-1.2 GeV for the correspondence to the simulation 
of the number of events. 




Fig. 2. Energy dependence of A and -P^. Left figure shows the magnitude of A. In the right figure, the 
two lines show the oscillation probabilities with 8 = 0° and 180°, respectively. 

In fig.2(left), one can see that the maximal value of A reaches almost about 0.2 and the probabilities 
for 6 = 0° and 180° have a difference about 0.4 around E ~ 0.2 GeV. In addition, the coefficient A is 
almost positive over the region E > 0.5 GeV and we expect a large CP dependence in total events 
accumulated in this region. If the coefficient A oscillates between positive and negative values, they 
contribute destructively and the CP phase effect becomes small in the total events. Therefore, the fact of 
A > in whole region of E > 0.5 GeV is very useful for the measurement of the CP phase effect. This is 
one of the important points of this paper. 

In fig. 2 (right), we can see that the probabilities in the case of 5 = 0° and S = 180° have a difference 
about 0.3 in the energy region E = 0.4-0.6 GeV. This can be understood by the magnitude of A for the 
corresponding energy region in fig. 2 (left). Up to now, the CP dependence in — > oscillations is 
neglected in a lot of works. Probably, one reason for this is that depends only on cos 5 and the CP 
dependence disappears in the case of 8 — 90° and 270°. Another reason is that the analysis satisfying the 
conditions E < 2 GeV and L > 2000 km is limited. In this letter, we found that there is a possibility for 
measuring 6 by using the channel of — ► oscillation under the above conditions. 

In order to understand the features in fig. 2, let us consider the approximate formula of A. The expression 
for ||HJ| is given by 

4A21A31S12C12S13C13 . A e L . A h L A32-L . . 

A — 7 — 7 sin - sin — - — cos — - — , (11) 

A* A f,. 2 2 2' v ' 
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where A h = Am 2 h /{2E), A e = Amj/(2E) and Ami (Am 2 ) is mass squared difference in matter 
The concrete expression for is given by 



(A 31 cos 29 13 - a) 2 + A 31 sin 2 28 13 ~ A 31 - a. (12) 

The second equality approximately holds in the case of sin #13 <C 1. We also obtain Ai by the replacements, 
A h -> An, A 3 i -> A 2 i and 9 13 -> 9 12 as 



(A21 cos 29x2 -a) 2 + A 2 i sin 2 29 12 ~ a. (13) 

The second equality approximately holds in the case of a ^ A21 = Am 21 /2E. If we substitute p = 
3.3g/cm 3 and the electron fraction in the mantle Y e — 0.494, we obtain a = 1.2 x 10~ 4 and Ag ~ a 
becomes a good approximation for E ^> Am 21 /2a = 0.3 GeV. Furthermore, at the baseline L = 3000 km, 
we obtain 

A/L aL 

sin ~ sin — ~ 0.8. (14) 
Therefore, the last two terms included in (jTTJ) are rewritten as 

. A h L A 32 L . (A 31 -a)L A 31 L 

2 sm cos ~ 2 sm cos 

2 2 2 2 

= |sin^A 31 i-^ -sin^j, (15) 

and they oscillate within the limits given by 

. aL A h L A 32 L . aL 

— 1 — sm — < 2 sm cos < 1 — sm — . (16 

2 2 2 2 K J 

Note that this takes value from —2 to 0, namely always negative, if sin ^ ~ 1. As a result, we can roughly 
estimate A as 

< < 4A 21 A 31 04^ ^ 013 

~ ~ a(A 3i - a) a E K ' 

in £ > 0.5 GeV and can understand the reason for being always positive in fig. 2 (left). We can also explain 
that the coefficient A is inverse proportional to the energy and decreases according to the increase of the 
energy. Furthermore, we find the position of the peak of A as 2 • 1.27(Am§ 1 — aE)L/E = (2n + 3/2)tt 
(n = 0, 1, 2, • ■ •) from fI5|l. If we substitute Am§ x = 2.2 • 10~ 3 eV 2 and L = 3000km, we obtain the peak 
energy as E = 1.27Am| 1 L/(2n + 7/4)7r = 5.3/(2n + 7/4). The peak with highest energy is obtained at 
E = 3 GeV by substituting n — 0. Thus, A has large value and is almost positive over the energy region 
about E = 0.5-5 GeV. Here, we also note the 9\ 3 dependence of A. As seen from (fTTj) . the CP phase effect 
becomes small according to the decrease of #13. 

Finally, let us briefly comment on the value of Jcp for — > v e oscillations as preparation of combined 
analysis using both and 24 events in the next section. We have investigated how Tcp for — > v e 
oscillations changes according to the values of E and L in the previous paper [12]. See fig.l (left) in 
ref.[12]. As a result, we found that the value of Icp takes nearly maximal value also in the region E < 2 
GeV and L > 2000 km for — + v e oscillations. This means that the absolute values of S' Te and S'^ e have 
an same order, namely the solar term is as large as the atmospheric term in this region. In the large Icp 
region, we can expect that the 5 dependence of the number of events also becomes large. 
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3 Estimation of Signal from Leptonic CP Phase 



Next, let us consider an experimental setup and estimate the CP phase effect. In order to utilize the S 
dependence of events in addition to that of v e events, we would like to choose the experimental setup, 
which satisfies the condition E ~ 0.5-2 GeV and L > 2000km obtained in sec. 2. As such an example, 
we consider the 4MW beam and the lMt water Cherenkov detector, which are the same used in the 
JPARC-HyperKamiokande experiment [23], but we take L = 3000 km as the baseline length. In fig. 3, 
we numerically calculate the signal of disappearance (left) and v e appearance (right), namely the total 
number of events distinct from the background noise, obtained in the above experimental setup within 
the energy range E = 0.4-1.2 when the CP phase 8 changes from 0° to 360°. We assume here only the 
neutrino beam data acquisition for two years and we use the same parameters as in fig. 2. We also give 
the statistical error within the 2-a level in fig. 3. We use the globes software to perform the numerical 
calculation 12311251. 
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Fig. 3. CP dependence of disappearance and v e appearance signal. Left and right panels show the 
and v e signal. The statistical error is also shown within the 2-<r level. 

Fig. 3 (left) shows that the number of z/ M events becomes maximal around <5 = 0° and minimal around 
8 = 180° as expected from the oscillation probability in fig. 2 (right). We obtain the number of events 
between these values in the case of 8 = 90° or 270°. We find that the number of events changes about 
250 due to the CP phase effect. On the other hand, fig. 3 (right) shows that the number of v e events 
becomes maximal around 5 — 45° and minimal around 8 = 225°. We also find that the number of v e 
events changes about 300 due to the CP phase effect. 

Next, let us explain the merit of the combined analysis of the total events in these two channels compared 
with the analysis of v e events alone. We assume <5 = 90° as the true value of the CP phase. Then, about 120 
v e events are expected and we obtain the information on the value of 8 as 0° < 8 < 15° or 75° < <5 < 105° 
or 350° < 8 < 360°. If the standard model is correct and the unitarity holds, about 1260 events are 
expected. From this result, we obtain the information 60° < 8 < 130° or 235° < 8 < 300°. If we combine 
the above information obtained by the two different channels, the extent of 8 is limited to 75° < 8 < 105°. 
Thus, there remain three allowed regions in the analysis of v e events alone. On the other hand, one allowed 
region is chosen in the combined analysis of two channels. Furthermore, if the number of events are 
largely different from 1260, this is considered to be the signal of some new physics beyond the standard 
model. 

In fig. 4, we also show the simultaneous fit of S13 and 8 by assuming the true values, sin 2 26>i3 = 0.16 
and 8 — 90°. Other parameters are the same as the values used in fig. 3. We define Ax 2 c and Ax 2 c+l/fl as 



A*2. = 
Ay 2 r 



N, 



true\2 



J\j~true 



N* 



(18) 
(19) 



where N V( and N^ ue represent the number of Vi events calculated by using the test values and true values 
respectively. In the left panel, the contours for 1, 2 and 3-cr C.L. with A\l c = 2.3,6.18 and 11.83 are 
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drawn as solid, dotted and dashed lines. In the same way, the contours for Ax^ e _|_„ are drawn in the right 
panel. In our analysis, only the statistical error is included for the simplicity. Here, we also use the globes 
software to perform the numerical calculation [24] [23] . 
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Fig. 4. Simultaneous fit of S13 and S by assuming the true values, sin 2 29\^ = 0.16 and 8 = 90°. Left and 
right panels are drawn assuming the measurement of v e events alone and both and v e events. The 
true values are represented by x in these figures. The statistical error is also shown as solid, dotted and 
dashed contours for 1,2 and 3-cr (A^ 2 = 2.3,6.18 and 11.83) respectively. 

Comparing the left and right panels, we also found that some allowed regions are excluded by the additional 
measurement of events. Another notable point is that the value of A\ 2 seems to remain unchanged 
even in the case for small sin 2 26*13. We can understand this result by the following reason. The appearance 
probability for to v e transition does not change largely because 6*13 independent term (the solar term) 
is comparatively large in such low energy region and is not so affected by the value of #13. Therefore, the 
value of Ax 2 does not change largely by the decrease of #13. 

Here, let us emphasize again on the importance of the simultaneous utilization of the — > and 
— > v e oscillations. This provides the variable possibility for the measurement of the CP phase in 
future experiments. We need to investigate further, including the parameter ambiguity problem [HHHH], 
background, systematic error and spectral dependence [21]. We also comment on the BNL-HS experiment 
[2"5] . In this experiment, L — 2540 km is taken as the baseline length and the peak energy of the neutrino 
beam is planned about E = 1-2 GeV. This setup satisfy the conditions E < 2 GeV and L > 2000km 
derived in this letter and is similar to our test experimental setup, so we expect the observation of the CP 
phase effect also in the BNL-HS experiment in the case of large 6*13. 

4 Summary and Discussion 

In summary, we have investigated the possibility for the measurement of the CP phase 8 by using — ► 
oscillations. The measurement of 8 in this channel have been considered to be difficult because the 
probability P w depends only on the cos 8 term and the CP dependence disappears in the case of 8 = 90° 
and 270°. We have used Icp and |A| (numerator of Iqp) for this channel and have found the CP phase 
effect becomes large in the region L > 2000km and E < 2 GeV beyond our expectation. As a result, we 
have shown that the difference between the probabilities for 8 — 0° and 180° reaches about 0.4 in the case of 
large 6*13. In addition, A is almost positive definite over the region E > 0.5 GeV and we expect a large CP 
dependence in total events accumulated in this region. As an example, we have used an experimental 
setup and have demonstrated that the allowed region is limited to one by combined analysis of v e and 
events although there remain three allowed regions by the analysis of v e events alone. Furthermore, this 
channel opens the window for a verification of the unitarity in three generation and exploration of new 
physics beyond the standard model. 
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